Dressed bound states for attosecond dynamics in strong laser fields 
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We propose a theoretical approach for the interpretation of pump-probe measurements where an 
attosecond pulse is absorbed in the presence of an intense laser pulse. This approach is based on 
abstractly defined dressed bound states, which capture the essential aspects of the interaction with 
the laser pulse and facilitate a perturbative description of transitions induced by the attosecond 
pulse. Necessary properties of dressed bound states are defined and various choices are discussed 
and compared to accurate numerical solutions of the time-dependent Schrodinger equation. 



INTRODUCTION 

The ionization of atoms and molecules in the field 
of an intense laser pulse is one of the most fundamen- 
tal phenomena in the interaction of light with matter. 
Despite decades of intensive experimental and theoret- 
ical research, there are still open conceptual questions, 
highlighted by recent experiments [H, 0|- A key ques- 
tion is whether one can define the instantaneous ioniza- 
tion probability in the presence of an ionizing field. A 
closely related question is to which extent the detach- 
ment of an electron by multiphoton absorption or quan- 
tum tunneling can be described by an ionization rate 
(see, for example, @, 0|)- Such a description is com- 
mon, but it can be made rigorous only in stationary 
or quasi-stationary situations. Two such cases are the 
static and time-periodic external electric fields. In static 
fields, the system decays exponentially at a well-defined 
rate. In time-periodic fields, Floquet theory allows defin- 
ing cycle-averaged decay rates (see e.g. [H or @ for a 
recent application). In both situations, real and imagi- 
nary parts of complex energies give energy positions and 
decay widths. Non-Hilbert space resonance states can 
be associated with the complex energies, either by us- 
ing complex scaling 0, [1[ or by iteratively adapting the 
boundary conditions (Siegert boundary conditions) [f|. 
None of these approaches, however, provides a fully sat- 
isfactory framework for investigating the interaction of 
few- to single-cycle ionizing pulses with matter. At the 
same time, the establishment of attosecond absorption 
spectroscopy 0, [1(3] and the recent advances in the gen- 
eration of optical waveforms [H, [H| call for a deeper un- 
derstanding of dynamics unfolding within a laser cycle. 
Such understanding implies the development of approxi- 
mate analytical theories, and the aim of this article is to 
make a step towards creating a useful theoretical frame- 
work. 

Our work is closely related to the development of adi- 
abatic or quasistatic approximations to strong-field ion- 
ization, such as described in [12| and the papers cited 
therein. In contrast to that research, we make no attempt 
to analytically describe strong-field ionization. Instead, 



we search for an analytical description of resonant tran- 
sitions that an attosecond pulse of ultraviolet (UV) or 
extreme-ultraviolet (XUV) radiation drives in the pres- 
ence of a strong near-infrared field. 



DRESSED BOUND STATES 

Let us consider a quantum system (atom, ion, molecule 
etc.) described by a Hamiltonian Hq. In the absence 
of any external field, this system has a set of orthonor- 
mal bound stationary states that satisfy the stationary 
Schrodinger equation: 



H |n) = ei 0) |n) 



(1) 



A typical attosecond measurement consists in letting the 
system interact with the electric fields of an intense laser 
pulse El(£), which may be regarded as a pump pulse, 
and a relatively weak probe pulse (or a train of pulses) 
Eprobe(i)- An XUV probe pulse may be as short as a few 
tens of attoseconds 13|, [l4j . Varying the delay between 
the two pulses with an attosecond accuracy and observ- 
ing the outcomes of their interaction with a quantum 
system provides valuable insights into phenomena trig- 
gered and steered by these two pulses [15| . If a significant 
fraction of atoms or molecules is ionized within a single 
half-cycle of the laser pulse, the interaction with the pulse 
is highly non-perturbative, so that an accurate descrip- 
tion of this interaction calls for a numerical solution of 
the time-dependent Schrodinger equation (TDSE). We 
assume that such a solution \tjj-^(t)) is available, and one 
of our goals is to interpret strong-field dynamics in terms 
of some dressed states. These states, if introduced prop- 
erly, may serve as a convenient starting point for building 
approximate models describing the interaction with the 
probe pulse. Before we discuss different approaches to 
defining such dressed states, let us list a few requirements 
that they should fulfill. 

Let \<f n (t)) be a bound dressed state associated with 
a stationary state |n). This means that \tp n (to)} = |ri) 
at some moment to before the interaction with external 
fields. As a first requirement, we demand that dressed 
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states should form an orthonormal set of wave vectors at 
all times: 



(<Pm(t)\ip n (t)} = 5 n 



(2) 



Orthonormality is not only convenient for computations, 
but it is also needed for the physical interpretation of 
the states: only for orthonormal states does finding a 
system instate \ip n (t)) with probability \{tp n (t)\' l pL(t))\ = 
1 imply that all other probabilities are 0. 

A strong laser field can significantly ionize a quan- 
tum system. As a free electron cannot absorb a photon, 
the probe pulse predominantly interacts with bound elec- 
trons. Therefore, it would be desirable to have dressed 
states that serve as a basis for describing the bound part 
of an electron wavefunction even in the presence of a 
strong external electric field. When this property is satis- 
fied, strong-field ionization should deplete bound dressed 
states without changing their shape in a fashion similar 
to the decay of quasi-stationary states. 

Except for losses due to ionization, each dressed state 
should fully account for distortions caused by the laser 
field. Therefore, we assume that the pump pulse induces 
no transitions between dressed states. This is the most 
important requirement that we place on dressed states, 

and it means that (^p m {t2) Uh(t2,ti) <Pn(ti)) must be 

negligibly small if m ^ n, where U^(t2,ti) is the propa- 
gator describing the interaction with the laser pulse (the 
exact definition of the propagator is given below). In 
particular, we expect that the interaction with the laser 
pulse brings the quantum system from its initial state 
\n) to the corresponding dressed state \<p n (t)) multiplied 
with a probability amplitude a n (t), the modulus of which 
accounts for the depletion by the laser pulse: 



a n (t) = Unit) U h {t,t ) 



while 



if m ^ n. 



(3) 



(4) 



As long as the above assumptions hold, the effect of the 
laser field on any dressed bound state can be expressed 
via the probability amplitudes: 



a n (h 



(5) 



and the time evolution of the system can be approxi- 
mated as 



#L(t2,il)«^|¥>n(t2)) 



a n (h) 



(<Pn(h)\ 



(6) 



Note that the approximate time evolution is no longer 
unitary as it accounts for ionization by losses from the 
dressed-state space. 



Assuming that time-dependent wave vectors \ip n (t)) 
satisfying all the above requirements exist, let us use 
them to develop a general perturbative description of 
the interaction of a laser-dressed system with a weak 
high-frequency probe pulse. For simplicity, we assume 
that the pump (laser) and probe (UV or XUV) pulses 
are polarized along the z-axis. In the dipole approx- 
imation, the length-gauge Hamiltonian takes the form 

H(t) = H + %{t) + V pTO be(t) With 



Mt) = eE L (t)Z, 
Vp TO be(t) = eE pro b c (t)Z , 



(7) 
(8) 



where e > is the electron charge and Z is the dipole op- 
erator. The interaction with the laser pulse is described 
by 



ih- |V>l(*)> = (H (t) + \Mt)) (9) 
or, in the operator form, 

ih^U h (t,t ) = (Ho + V L (t))U L (t,t Q ) (10) 

with IV'l(O) = Ujjitjto) \i) for an initial state A stan- 
dard approach to developing a perturbation theory with 
respect to the probe pulse consists in writing the com- 
plete TDSE 



ihj t \m) = H(t) \m) 



(ii) 



in the following integral form: 



U(t,t ) - U L {t,to) 



dt'U(t,t')V pIobe (t')U h (t',t ). 

(12) 

Here, U(t,to) is the propagator associated with the full 
Hamiltonian H(t): 



ih^-U(t,t )=H(t)U(t,t ) 
at 



(13) 



and IV^)) = U(t,to) \i) for the initial condition |^(t )) = 
\i). Furthermore, if only single-photon processes play a 
role in the interaction with the probe pulse, then the un- 
known operator U(t, t') on the right-hand side of Eq. (fT2|) 
may be approximated with Ui,(t,t'), which we consider 
to be known. This yields 



|V>(*)> « \Mt)) - r / dt'U h (t,t')V piohe (t')\Mt'))- 

(14) 

The probe pulse may cause bound-bound, as well as 
bound-continuum transitions. In this paper, we focus 
on transitions between bound states, leaving the direct 
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photoionization by the probe pulse aside. In this case, 
Eq. (fT4)) can be written as 

\m) « \Mt)) - j E f dt ' (wo mo) x 



fl(t') %Tobe(t') <p m (t'))(<p m (t')\llJ L {t'))}. (15) 



Let us evaluate the probability amplitude that a sys- 
tem, initially prepared in state \ip(to)) = \i), will be found 
in state \ip n ^i(t)) at a later time t. Multiplying both sides 
of Eq. (fT5|) with (ip n (t)\ from the left, keeping in mind 
the fact that [iphit)) s=s ai(t) \ifi(t)), and using Eq. ([5]), 
we obtain 

a m (t) = (ip n {t) \u{t,t )\ = (Vn(*)|^(<)> « 

-~ea n (t) [ dt'E plohe (t')^Z ni (t') (16) 



for n^i, where 



(17) 



is the transition matrix element between dressed states. 
Eq. (fT6"|) is the desired expression, where the interaction 
with the probe pulse is accounted for in first-order per- 
turbation theory, while the effect of the strong laser pulse 
is represented by quantities related to dressed states. 
Based on Eq. (TT6"|) . we can evaluate quantities observable 
in measurements. For example, the probability that, at 
some final time if after the interaction, the system will 
be found in a bound stationary state |/} is given by 



p/ { =iw(tf))r 



^Uni(U) {f\<fn(tf)) 



(18) 



As an example that is more relevant to transient ab- 
sorption spectroscopy, the dipole response due to bound- 
bound transitions induced by the probe pulse can be eval- 
uated as 

d(t) = fat) |z| V(t)) ^^a* mj {t)Z mn {t)a ni {t). (19) 



So far, the dressed bound states were treated as an 
abstract mathematical concept. It is worth mentioning 
that without the requirement that dressed states should 
describe bound electrons, the wavevectors \ip n (t)) = 
Uh(t,to) \n) would exactly satisfy the rest of our require- 
ments, but such dressed states would not serve their 
main purpose, which is to facilitate the interpretation 
of bound electron dynamics in the presence of an intense 
laser pulse. We do not know if perfect dressed states sat- 
isfying all our requirements including their bound char- 
acter exist, but we can consider several approximations 



to dressed bound states and use our general theory to 
investigate their performance in comparison with exact 
solutions of the TDSE. We will consider four different 
definitions of dressed bound states. 

As a very crude approximation, one can attempt to use 
the unperturbed stationary states: 



(20) 



Another option, frequently encountered in literature, is 
to use time-dependent eigenstates of the Hamiltonian 
H(t) evaluated in the subspace of unperturbed bound 
states and normalized to satisfy Eq. ([2]): 



E 



f ( n\t) 



= e {a) 



(*) 



(21) 

In the context of strong-field ionization, these states are 
known under many different names, such as adiabatic 
states , quasistatic states 1171 . Il8| , phase-adiabatic 

states [19(, field-adapted states [20|, and adiabatic field- 
following dressed states 2l|, |22j . We will refer to them 
as "adiabatic states" . 

Yet another reasonable approach is to obtain dressed 
bound states by solving the TDSE in the basis of unper- 
turbed bound states: 



in— 

dt 



^ d) (*))=Ei 



with the initial condition 



H + V L (t) 



(V 



(22) 



|n). With this 



choice, we assume that virtual excitations into field-free 
continuum states do not influence the dynamics. How- 



ever, contrary to 



fn (t) ) , these states do not follow 

the laser field adiabatically, but rather allow for non- 
instantaneous response to the laser field. We will call 
them "dynamically dressed states" . 

The above three alternatives allow one to evaluate the 
dressed bound states without solving the TDSE in the 
full Hilbert space, although such solutions \iph(t)) are re- 
quired later to obtain the probability amplitudes a n (t) 
from Eq. ([3]). The last option that we are going to con- 
sider is different: we will evaluate dressed bound states 
by projecting solutions of the TDSE obtained in the 
whole Hilbert space ([9]) onto the subspace of unperturbed 
bound stationary states and orthonormalizing them with 
the aid of the Gram-Schmidt procedure: 



^\t)) = n t 1 £w) e i™> ( m 

Kn m 



(23) 



^Ht)) = \^\t)) /yj($\W<*\t)), (24) 



where the projector Pi(t) is defined as 

p l (t)= <p<r\t))/ v <r\t) 



(25) 
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We will call 



Mi 



fn \t)) "projected dynamical states", as 
they result from calculations that account for all electron 
dynamics including strong-field ionization, but the out- 
comes of these calculations are projected onto the field- 
free bound-state content. Note that, evaluating the pro- 
jected dynamical states, the projected dressed ground 
state is only normalized without any orthogonalization, 
the first excited dressed state is forced to be orthogonal 
to the dressed ground state and so on. 



NUMERICAL RESULTS 

Unless stated otherwise, we use atomic units (at. u.) 
in this section: h = e = m e = 1, where the units of en- 
ergy and length are 1 Hartree = 27.21 eV and the Bohr 
radius (5.29 x 1CP 11 m), respectively. One atomic unit 
of the electric field is 5.412 x 10 11 V/m. To compare the 
different approaches to defining dressed bound states and 
to illustrate the power of our analytical theory, we solve 
the TDSE numerically for one electron in one spatial di- 
mension with a soft-core model potential: 



V(z) 



1 



(26) 



sc 



With the soft-core parameter being equal to asc = 
0.3, the deepest two energy levels have energies e$ = 
-1.75 = -47.49 eV and e^ 0) = -0.41 = -11.11 eV. The 
TDSE was solved on a grid with a step of Az = 0.1 
atomic units in a box as large as 819.2 atomic units. We 
used the first five bound states to investigate various def- 
initions of dressed states. The most excited of these un- 
perturbed states has a binding energy of £4°' = —0.077 = 
-2.10 eV. 

We define the light pulses via their vector potentials: 



A(t) 



-Anax cos 2 (t/T) sin(wt) for \t\ < ttT/2 and 



A(t) = for \t\ > ttT/2. Here, A max = S max /w is the 
peak value reached by the envelope of the vector poten- 
tial, -E max is the peak value reached by the envelope of 
the electric field, u is the central frequency, and T is 
related to the full width at half maximum (FWHM) of 
the pulse intensity as T = FWHM/ [2 arccos (2~ 1/4 )] . 
Given a vector potential A(t), the electric field is evalu- 
ated as F(t) = —A'(t). In our simulations, we used a 3.5- 
femtosecond laser pulse (T L = 126.78 = 3 fs) and a 300- 
attosecond XUV probe pulse (T probe = 10.84 = 0.26 fs). 
The central wavelength of the laser pulse was set to 760 
nanometers (wl = 0.06 = 2.48 fs -1 ). The central fre- 
quency of the XUV pulse was chosen to be resonant with 
the transition between the first two stationary states: 
Wprobc = ^ 0) - 4 0) = 1.34 = 36.38 eV. Although we 
make no attempt to model a realistic atom, the excita- 
tion from |0) to |1) may be viewed as analogous to a 
dipole-allowed transition from an inner-shell orbital of a 
multi-electron atom to an unoccupied orbital. 




0.02 



FIG. 1. The probabilities of being in variously defined dressed 
states \<pi(t)}, provided that the quantum system is initially 
prepared in state |1). The modulus of the laser field is 
shown with black dots. The two panels correspond to dif- 
ferent peak intensities I ma x of the laser pulse: a) E msiK — 
0.02 atomic units (/ max = 1.4 x 10 13 W/cm 2 ) and b) £ max = 
0.04 atomic units (J max = 5.6 x 10 13 W/cm 2 ). 



In Fig. [TJ we compare the probabilities |a„(i)| 2 = 
I (<p n (i)\i/ji,(t)) j 2 evaluated using the four different def- 
initions of dressed bound states introduced in the pre- 
vious section. In these simulations, we set n = 1, 
that is, each simulation begins with the electron be- 
ing in first excited state |1). For the peak electric field 
B max = 0.02 = 10 10 V/m, the ionization probability is 
as little as 0.3%. If we knew the perfect dressed bound 
states, then, for such a small ionization probability, our 
model system would largely remain in \ipi(t)}, so that 
|ai(i)| 2 would only take values between |eti(to)| 2 = 1 and 
|ai(if)| 2 = 0.997, decreasing more or less monotonously 
as the initial state is slightly depleted. Fig. [1] illustrates 
the well-known fact that the unperturbed bound states 
are very far from this ideal: |di (i)| 2 reaches values be- 
low 0.95. This happens mainly because the initial state 
is distorted by the external field, which reduces the value 
of I (n\ipi^(t)) | 2 . So, even though |a^"'(i)| 2 can be inter- 
preted as the probability of being in the unperturbed sta- 
tionary state, it should not be given the physical meaning 
of being in some distorted but bound state that, after the 
interaction, turns into the initial unperturbed stationary 
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defined by 



state. 

In this respect, dressed state 
Eq. (PHI) , performs much better, but it is not perfect 
either. Even though tp^\t)^j accounts for the distor- 
tion of the initial state in the external field, |aj (i)| 2 
in Fig. [l}i still has pronounced minima at times where 
E\,{t) has its zero crossings, and |a^(£)| 2 reaches values 
that are significantly smaller than the final probability 
of being in state The main reason for this is that 
| fi^ CO ^) is a stationary solution, while the average ve- 
locity of a bound electron can become relatively big in 
an intense laser pulse. As this velocity has maxima at 



<P l i\t)\Mt) 



has 



zero-crossings of the electric field, 

minima at these times. This can also be interpreted as a 
breakdown of the adiabatic (quasistatic) approximation 
(23j . Also the small delays between the extrema of (t) 
and the minima of |a^(t)| 2 , visible in Fig. [1^, point out 
some non-adiabaticity of the electron response. 

The intuitively expected step-like decrease of ai(t)| 2 , 
as the laser field depopulates the initial state, is best 



reproduced by the dynamically dressed states 



and the projected dynamical states 
by Eqs. {22} and {24 
and 



^0) 
defined 



respectively. However, unlike 
ifn '(t)V these dressed states do not nec- 
essarily turn into \n) once the laser pulse is gone. In 
other words, \a[ d ' p \t{)\ 2 ^ | (l|V>L(*f)) | 2 at a final time 
t{. As long as the discrepancy is small, it is legitimate to 
interpret |ai(i)| 2 as the probability of being in the first 
excited state dressed by the laser field. In these simula- 
tions, this discrepancy is much smaller for |a^(i)| 2 , as 
compared to |a^(t)| 2 , which is especially easy to see in 
Fig. [Tb. 

The mere fact that, for a certain choice of dressed 
states, |a n (i)| 2 does not show the intuitively expected 
step-like behavior, does not necessarily mean that these 
dressed states should not be used in our analytical result 
(IT6T) for the interaction with a probe pulse. In Fig. [21 
we plot the results of simulations where our model atom, 
initially prepared in state |1), interacts with the near- 
infrared laser pulse 2?l (0 and a delayed attosecond XUV 
pulse E pTO bc(t — Tprobo). We plot the probability to find 
the atom in state |0) at the end of the simulation, nor- 
malized to the probability of the same transition in the 
absence of the laser pulse. Applying our model, we found 
that care should be taken implementing Eq. (fT5|) . If a 
highly excited bound state is strongly depleted by the 
laser pulse, a n {t') in the denominator of Eq. ([T6|) may 
become very small. As a result, even small inaccuracies 
in the approximations that we made deriving Eq. (|16[) 
may result in big errors. To avoid this, we restrict the 
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FIG. 2. The effect of the laser field on the probability to 
find our model quantum system in state |0) at a final time tf, 
the initial state being |1). The XUV pulse, which is mainly 
responsible for the transition, arrives at a time T pr obc. The 
probabilities are normalized to the transition probability in 
the absence of the laser field. The approximate probabilities 

Pol j P0X1 P01 j an d P01 are evaluated with Eq. (|27[) . using 
different definitions of dressed bound states. The exact tran- 
sition probabilities Poi DSE , obtained by numerically solving 
the TDSE, are shown with filled circles. The peak intensities 
of the laser pulse are the same as those in Fig. [T] -B max = 0.02 
for panel a) and -E max = 0.04 for panel b). 



number of states used in Eq. (|T8 



Pfi 



max{i,/} 

a ni(t{) (f\<Pn(ti)) 



(27) 



The filled circles in Figs. [2^, and[UD represent the results 
of solving the TDSE on a grid. Comparing these "ex- 
act" transition probabilities with the predictions of our 
analytical model, we clearly see that projected dressed 



states 



Mi 



tfn '(t)), defined by Eq. , give the most accu- 
rate results. At the same time, the unperturbed states 
<Pn^\ = \n), which are our most primitive dressed states, 
also perform surprisingly well. For the calculations with 
the higher intensity, presented in Fig. [2p, these states 



even outperform both the adiabatic states 



and 



the dynamically dressed states 



(t) ) in terms of the 
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agreement with the TDSE calculations. For this high in- 
tensity, the biggest discrepancy between the exact and 
analytical results is found for the dynamically dressed 



states 
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FIG. 3. Fhe effect of the laser field on the transition proba- 
bility from state |0) to state |1). For details, see the caption 
of Fig. [5] and the text. 

In Fig. [21 we present similar calculations for a case 
more relevant to possible experiments — in this case, the 
XUV pulse excites our model atom from its ground state. 
The figure depicts the probability to find the atom in 
state |1) after the interaction with the XUV and laser 
pulses. In contrast to Fig. [21 the unperturbed states 
Pi^it)/ fail to describe the modulation of the transi- 



tion probability by the laser field. Both 



ip[ a) (t)) and 



fi W ) yield transition probabilities that agree with the 
TDSE calculations. The corresponding discrepancies are 
comparable for the more intense laser pulse (Fig. [3jj) , but 

using ^^(i)^ clearly gives more accurate results if the 
depletion of state |1) is weak (Fig. fS^)- Thus, in most 
cases, the projected dynamical states ifin \t)^, defined 

by Eq. (fMf , best satisfy the requirements that we placed 
on dressed bound states. 

Having identified the most promising definition of 
dressed bound states, we can interpret the effect of the 
laser pulse on the probabilities of XUV-induced transi- 
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FIG. 4. The squared modulus of the dipole transition ma- 
trix element between dressed states jv^cT'W^) an( i 
as defined by Eq. (|17[) . j5 max is the peak value of the elec- 
tric field of the laser pulse. The modulus of the laser field is 
schematically shown with black dots. 



tions. The overall decrease of poi with r pro b G in Fig. EJd, 
as well as the increase of pio in Fig. [3t> are clearly due 
to the depletion of the first excited state by the laser 
field. The periodic modulations present in all the fig- 
ures deserve more consideration. In principle, they may 
be a consequence of two effects: the Stark shift, which 
detunes the transition frequency away from the central 
frequency of the XUV pulse, and the modulation of the 
dressed transition matrix element (|17[) . In the present 
example, we find that the dominant contribution is due 
to the effect of the laser field on the dressed matrix ele- 
ments. To see this, we plot |Zio(i)| 2 , evaluated with the 



aid of 



>b 



fn (t)h f° r different intensities of the laser field: 
Fig. |4j For the central half-cycles of the laser pulse, the 
modulation depths of |Zio(i)| 2 are very close to those of 
transition probabilities in Figs. [2] and G3 A strong field 
tends to reduce |i?io(i)| 2 because it displaces the electron 
wavefunction in the first excited state, while the ground 
state in the present example has a very small polariz- 
ability. As a consequence, the overlap between the two 
dressed states decreases, which results in a smaller value 
of the dipole transition matrix element between them. 

The relatively fast oscillations of |Zio(i)| 2 at the tail 
of the laser pulse are mainly due to the multiphoton ex- 
citation of state 1 3) from state |1). The energy difference 



between these two states is e 



(0) M 



-0.113 + 0.408; 



0.29 atomic units; since the transitions from both these 
states to |0) are dipole- allowed, this leads to quantum 
beats with a period of 27r/0.29 = 21 at. u. = 0.5 fs. 
These quantum beats do not, however, appear in the 
probabilities of the XUV-induced transitions presented 
in Figs. [2] and [3] because the XUV pulse, which is reso- 
nant with the transition between states |0) and |1), does 
not have enough bandwidth to drive transitions between 
states |0) and |3). This fact is successfully accounted for 
by our analytical theory. 
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CONCLUSIONS AND OUTLOOK 

We have theoretically investigated how an ionizing 
few-cycle near-infrared laser pulse affects single-photon 
bound-bound transitions driven by an attosecond pulse 
of extreme ultraviolet radiation. Our approach is based 
on the assumption that even when the laser field is strong 
enough to significantly ionize a quantum system, one can 
describe bound electrons with some dressed bound states, 
which we introduce in an abstract way by a set of simple 
requirements listed at the beginning of section „ With 
such dressed states, we have obtained relatively simple 
analytical results for the perturbative interaction with a 
probe pulse, such as Eq. For a particular model 

problem, we have systematically compared four different 
kinds of dressed bound states and found that the most 
accurate results are generally obtained with the projected 
dynamical states ([24|l . These states are easily evaluated 
from numerical simulations of the interaction with the 
laser pulse, and they make no assumptions about the 
adiabaticity of the response to the laser pulse. 

As the main observable in our numerical examples, we 
have chosen the probability to find a single-electron atom 
in a certain state after the interaction with light pulses. 
This was done to test our analytical theory in a possibly 
simple and convincing way. Experimentally, it is easier to 
measure the transient absorption of an attosecond pulse 
in the presence of an ionizing laser field [l(|. It should 
be straightforward to employ dressed bound states to 
theoretically investigate attosecond transient absorption 
measurements. In this paper, we have only made a first 
step in this direction by writing an explicit expression 
for the dipole response associated with an XUV-driven 
laser-dressed bound-bound transition: Eq. (|19j) . 

The main purpose of this work has been to define gen- 
eral properties for dressed states in strong laser-atom in- 
teractions and to explore their existence. We have shown 
that quantum states with properties close to those of ide- 
alized dressed bound states can be found. From the fact 
that, overall, the projected dynamical states ip^ perform 
best, while all other states show severe shortcomings in 
various situations, we can draw the following preliminary 
conclusions: Not surprisingly, beyond the lowest intensi- 
ties, the unperturbed states are not suitable for the 
analysis of laser induced dynamics. Adiabatic states ip( a ) 
provide a good qualitative picture, but they fail when the 
velocity acquired by a bound electron in a strong field 
becomes significant, which is the case for systems with a 
reasonable polarizability. The dynamically dressed states 
showed us that admitting dynamics, but restricting 
it to the field-free bound states, may worsen rather than 
improve the agreement with the full TDSE: it appears 
that the restriction of the dynamics actually introduces 
serious artifacts, even when comparing to the simpler 
adiabatic states computed in the same subspace of field- 



free bound states. Finally, from the projected dynamical 
states ip( p ) we learn that the field-free continuous part of 
the spectrum influences the dynamics, but does not play 
a crucial role for the computation of transition matrix 
elements between bound states. These four examples are 
clearly not the only possible choices. For example, rather 
than restricting to a few field-free bound states, one may 
use adiabatic or dynamic states confined to a box, which 
would admit virtual continuum states without the need 
for a full solution of the TDSE, as in the projected dy- 
namical states. 

We have thus shown that also in situations where sig- 
nificant and strongly time-dependent ionization occurs, 
dressed bound states allow for a simple analysis of ex- 
perimentally observable dynamics and carry our picture 
of state-depletion and perturbative transitions between 
the dressed states into extreme time scales and intensity 
regimes. This finding is conceptually important. Prag- 
matic considerations how to further improve the dressed 
states and possibly find a compromise between accuracy 
and computational simplicity are left to future work. 
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